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Abstract 

String inspired models with millimeter scale extra dimensions provide a natu- 
ral way to understand an ultralight sterile neutrino needed for a simultaneous 
explanation of the solar, atmospheric and LSND neutrino oscillation results. 
The sterile neutrino is the bulk neutrino (i/g) postulated to exist in these 
models, and it becomes ultralight in theories that prevent the appearance of 
its direct mass terms. Its Kaluza-Klein (KK) states then add new oscillation 
channels for the electron neutrino emitted from the solar core. We show that 
successive MSW transitions of solar v e to the lower lying KK modes of vb in 
conjunction with vacuum oscillations between the v e and the zero mode of vb 
provide a new way to fit the solar neutrino data. Using just the average rates 
from the three types of solar experiments, we predict the Super-Kamiokande 
spectrum with 73% probability, but dips characteristic of the 0.06 mm extra 
dimension should be seen in the SNO spectrum. We discuss both intermedi- 
ate and low string scale models where the desired phenomenology can emerge 
naturally. 

PACS: 14.60.Pq; 14.60.St; ll,10.Kk; 



I. INTRODUCTION 

At present there appear to be three classes of experiments that provide evidence for 
neutrino oscillations: solar neutrino searches atmospheric neutrino data p[ and an ac- 
celerator search for oscillations by the LSND experiment ||. A simultaneous understanding 
of all these data seems to require the existence of an ultralight neutrino (beyond the three 
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known ones: v e , i/ T ), which must be sterile with respect to weak interactions. Within this 
four- neutrino scheme ||], the solar v e deficit is explained by v e — * v s (where v s is a sterile 
neutrino), the atmospheric v^/v e anomaly is attributed to — > u T , and the LSND || results 
are explained by the v e — oscillation predicted in the model. The heavier near-degenerate 
V/j, and v T are required by the LSND results to be in the eV range and can therefore share the 
role of hot dark matter. Exactly this same pattern of neutrino masses and mixings appears 
necessary to allow production of heavy elements (A ~ 100) by type II supernovae M. While 
qualitatively this neutrino scheme || seems to explain all existing neutrino phenomena, solar 
neutrino observations are now sufficiently constraining that the small-angle MSW v e — > v s 
explanation appears to be in some difficulty J7|, and seemingly one must go to some length 
|| in order to try to rescue this scheme. Although providing better fits to the solar data, 
even active-active transitions in a three-neutrino scheme do not give a quantitatively good 
explanation of all those data. The theoretical and phenomenological challenge then is to 
find a scheme which has an ultralight sterile neutrino and at the same time provides a fit to 
the solar neutrino data. 

Recently, motivated by string inspired brane models with large extra dimensions ||, we 
pointed out that there appears to be a way to achieve an excellent fit and rescue the 



apparently needed four- neutrino scheme by including a singlet, sterile neutrino in the bulk. 
The method provides nearly maximal u e vacuum oscillation with the lightest pair of KK 
modes, and also has MSW transitions to several of the other modes. We also pointed out 
that the model has its characteristic predictions for the v e survival probability and can be 
tested as new solar neutrino data on the neutrino energy distribution accumulates. It is the 
goal of this paper to elaborate on this proposal and discuss theoretical schemes that can 
lead to the desired parameters for the neutrinos. 

In section II we discuss the examples of brane-bulk models which lead to light neutrinos 
and the questions of naturalness of an ultralight sterile neutrino; in sections III and IV 
we consider two examples of models which have neutrino spectra desirable from the point 
of view of understanding neutrino data. In section V, we provide a fit to solar neutrino 
observations in the context of these models, taking matter effects into account exactly. We 
present our predictions for the recoil energy distribution for the Super-K data and annual 
variation of the flux. In section VI we summarize our conclusions. In Appendices A, B we 
provide some more details on the TeV scale as well as the local B-L models, and in Appendix 
C we comment further on the naturalness of ultralight sterile neutrinos in bulk-brane models. 



II. NEUTRINO MASSES IN MODELS WITH LARGE EXTRA DIMENSIONS 

One of the important predictions of string theories is the existence of more than three 
space dimensions. For a long time, it was believed that these extra dimensions are small and 
are therefore practically inconsequential as far as low energy physics is concerned. However, 
recent progress in the understanding of the nonperturbative aspects of string theories have 
opened up the possibility that some of these extra dimensions could be large |TT|,P| without 



contradicting observations. In particular, models where some of the extra dimensions have 
sizes as large as a millimeter and where the string scale is in the few TeV range have attracted 
a great deal of phenomenological attention in the past two years ||. The basic assumption 
of these models, inspired by the D-branes in string theories, is that the space-time has a 
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brane-bulk structure, where the brane is the familiar (3+1) dimensional space-time, with 
the standard model particles and forces residing in it, and the bulk consists of all space 
dimensions where gravity and other possible gauge singlet particles live. One could of 
course envision (3+d+l) dimensional D-branes where d-space dimensions have miniscule 
(< TeV^ 1 ) size. The main interest in these models has been due to the fact that the low 
string scale provides an opportunity to test them using existing collider facilities. 

A major challenge to these theories comes from the neutrino sector, the first problem 
being how one understands the small neutrino masses in a natural manner. The conventional 
seesaw [12 explanation which is believed to provide the most satisfactory way to understand 
this, requires that the new physics scale (or the scale of SU{2) R x U{1)b-l symmetry) be 
around 10 9 to 10 12 GeV or higher, depending on the Dirac masses of the neutrinos whose 
magnitudes are not known. If the highest scale of the theory is a TeV, clearly the seesaw 
mechanism does not work, so one must look for alternatives. The second problem is that 
if one considers only the standard model group in the brane, operators such as LHLH j 'M* 
could be induced by string theory in the low energy effective Lagrangian. For TeV scale 
strings this would obviously lead to unacceptable neutrino masses. 

One mechanism suggested in Ref. |D| is to postulate the existence of one or more gauge 
singlet neutrinos, z/g, in the bulk which couple to the lepton doublets in the brane. After 
electroweak symmetry breaking, this coupling can lead to neutrino Dirac masses, which are 
suppressed by the ratio M*/Mp£, where Mpi is the Planck mass and M* is the string scale. 
This is sufficient to explain small neutrino masses and owes its origin to the large bulk 
volume that suppresses the effective Yukawa couplings of the Kaluza-Klein (KK) modes of 
the bulk neutrino to the brane fields. In this class of models, naturalness of small neutrino 
mass requires that one must assume the existence of a global B-L symmetry in the theory, 
since that will exclude the undesirable higher dimensional operators from the theory. 

An alternative possibility |14| is to consider the brane theory to have an extended gauge 



symmetry which contains B-L symmetry as a subgroup. Phenomenological considerations, 
however, require that the local B — L scale and hence the string scale be of order of 10 9 GeV 
or so. The extra dimensions in these models could also be large. Indeed, it is interesting 
that if there were only one large extra dimension, R (and all small extra dimensions are 
~ M~ l ), the formula 

M 2 pe ~ MlR (1) 



leads to M* ~ 10 GeV if R~ mm [Q. While in these models, there is no strict need 
to introduce the bulk neutrinos to understand the small masses of known neutrinos, if we 
wanted to include the sterile neutrinos, one must add the ub- The high scale models may 
also have certain other advantages which we will see as we proceed. 

Regardless of which path one chooses for understanding small neutrino masses, a very 
desirable feature of these models is that if the size of extra dimensions is of order a millimeter, 
the KK excitations of the bulk neutrino have masses of order 10~ 3 eV, which is in the range 
needed for a unified understanding of oscillation data 0], as already noted. 
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III. TEV SCALE MODELS 



To discuss the mechanisms in a concrete setting, let us first focus on TeV scale mod- 
els. Here, one postulates a bulk neutrino, which is a singlet under the electroweak gauge 
group. Let us denote the bulk neutrino by vb{^iV)- The bulk neutrino is represented 
by a four-component spinor and can be split into two chiral Weyl 2-component spinors as 
u b = (x T y ~ The 2-component spinors x an d 4> can be decomposed in terms of 
4-dimensional Fourier components as follows: 

/ \ 1 1 ( n^V ■ ■ nny\ 

= + \ X+ ' nC ° S ~R~ + l X-,nSin— j . (2) 

There is a similar expression for 0. It has a five dimensional kinetic energy term and a 
coupling to the brane field L(x^). The full Lagrangian involving the vb is 

£ = WbI^Vb + kLHv B r(x, y = 0) + i J dy v B l(x, y)d^ BR {x, y) + h.c, (3) 

where H denotes the Higgs doublet, and k = hj^- is the suppressed Yukawa coupling. This 
leads to a Dirac mass for the neutrino |L3| given by: 



hv wk M* 

m = - (4) 

where v w k is the scale of SU(2)l breaking. In terms of the 2-component fields, the mass 
term coming from the fifth component of the kinetic energy connects the fields x+ with 0_ 
and X- with <f> + , whereas it is only the + (or v b ,r,+) which couples to the brane neutrino 
v e L . Thus as far as the standard model particles and forces go, the fields 0_ and x+ are 
totally decoupled, and we will not consider them here. The mass matrix that we will write 
below therefore connects only v eLl + n and X-,n- 

From Eq. |], we conclude that for M* ~ 10 TeV, this leads to m ~ 10~ 4 /i eV. It is 
encouraging that this number is in the right range to be of interest in the discussion of solar 
neutrino oscillation if the Yukawa coupling h is appropriately chosen. Furthermore, this 
neutrino is mixed with all the KK modes of the bulk neutrino, with a mixing mass ~ 
since the nth KK mode has a mass nR~ l = nfi, the mixing angle is given by \[2mRjn. 
Note that for R ~ 0.1mm, this mixing angle is of the right order to be important in MSW 
transitions of solar neutrinos. 

It is worth pointing out that this suppression of m is independent of the number and 
radius hierarchy of the extra dimensions, provided that our bulk neutrino propagates in the 
whole bulk. For simplicity, we will assume that there is only one extra dimension with radius 
of compactification as large as a millimeter, and the rest with much smaller compactification 
radii. The smaller dimensions will only contribute to the relationship between the Planck 
and the string scale, but their KK excitations will be very heavy and decouple from the 
neutrino spectrum. Thus, all the analysis can be done as in five dimensions. 

In order to make this model applicable to resolving the observed oscillation phenomena, 
we have to extend the model as has been noted in . Even if one wanted to understand the 
solar neutrino oscillation using this picture, one would have difficulty fitting all the rates in 
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Gallium, Chlorine and the water Cherenkov data while simultaneously explaining the recoil 
energy distribution in the Super-Kamiokande data. 

One way to do this would be to include new physics in the brane. We parameterize this 
in terms of an effective Majorana neutrino mass matrix in the brane: 



M 




(5) 



The origin of this pattern of brane neutrino masses will be discussed in Appendix A. In 
this section we concentrate on the effect of this matrix on the mixing of the bulk neutrinos 
with the brane ones. For this we will assume that m ^> 5^; as a result, the v^ T decouple 
and do not affect the mixing between the bulk neutrino modes and the u e , and in the 
subsequent analysis we consider the remaining modes. Their mass matrix in the basis 
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One can evaluate the eigenvalues and the eigenstates of this matrix. The former are the 
solutions of the transcendental equation: 
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The equation for eigenstates is 
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(7) 



(8) 



where we have used the notation ± for the left- and right-handed parts of the KK modes of 
the bulk neutrino in the two-component notation and dropped the L, R subscripts, the sum 
over k runs through the KK modes, and N n is the normalization factor given by 



A/-2 = i + m 2 n 2 R 2 + 



m c 



(9) 



From Eq. [7| we see that there are two eigenvalues m± near m when (5 ee C m <C yU, and these 
are given by m± ~ ^ ±m. These are the lowest two levels, and their mass difference square 
is given by ~ 2m5 ee eV 2 . From Eq. [8], we see that they are maximally mixed. Therefore, 
if S ee ~ 10~ 7 eV, then the transition between these levels can lead to vacuum oscillation 
("VO") of the solar neutrinos. This will be one of the ingredients of our new solution, as we 
see below. 
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IV. LOCAL B — L SYMMETRY MODELS: 



A second way to achieve the same phenomenology is possible using a much higher string 
scale. In this class of models |TJ| , one postulates that the theory in the brane is left-right 
symmetric so that it contains the B-L as a local symmetry. The gauge group of the model is 
SU{2) L x SU{2) R x U(1)b-l with field content for leptons given by left and right doublets 
ipL,R = iyA)L,R and in the Higgs sector bidoublet <fi = (2,2,0), doublets Xl,r- As in the 
previous case we choose a single bulk neutrino. We then impose the following Z2 symmetry 
on the Lagrangian under which xl ~* ~Xl, ^b,r — > —vb,r and all other fields are invariant. 
Note that under this symmetry, the 5th coordinate y — > —y. The invariant Lagrangian is 



/ _ 

+ 1pLH>R + T7L/2 bPeRXR + 1peLXL]vB(x,y = 0) 



M 



+i J dy v BL (x,y)d 5 p BR (x,y) +h.c, 



(10) 
(11) 



where / and hg are Yukawa couplings and Xl,r = i> T 2X*L r- We then break the right-handed 
symmetry with (x° R ) = v R , while at the same time keeping (xl) = 0. We expect v R to be of 
order of the string scale, M*. The Lagrangian involving the electron neutrino and the bulk 
neutrinos then becomes 

v 2 p 
C = Ij^vrVr + mv eL v eR + au eR u BL (x, y = 0) + i J dy v B l{%, y)d 5 is BR (x, y) + h.c, (12) 



where a ~ hl M*J l R an< l m = Mwfc- In this section we discuss only the v e — v B sector and 
address the full three generation mixing in Appendix B. In what follows, the Majorana mass 

2 

of the v R is denoted by M, i.e., M = fj*h- This leads to the mass matrix of Eq. [13] which 
mixes the brane neutrinos with the KK modes of v B . For (x° R ) = v R ^> (</»), /x, which we 
assume, u eR has a Majorana mass, M, which is much bigger than any other masses in the 
theory (ignoring real superheavy KK modes), and v eR decouples. One can then write an 
effective Lagrangian at low energies (i.e., E -C v R ) using the seesaw mechanism. The effective 
mass matrix for the light modes can be written down using the same notation for the KK 
modes of v B as in the previous section. In the basis [y e Li v B l +> u bl +> u br - > u bl +j v br -■>'■>') 
it is given by: 



M = M 
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(13) 



It is easy to see that the lowest eigenvalue of this matrix is zero. The transcendental equation 
describing the rest of the eigenvalues is 



m 2 a 2 nm n irm T 

1 -cot 

M M /j, fj. 



(14) 
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The next lowest eigenvalue solution of this equation is mi ~ m ^ a ■ For m, a ~ 1 — 5 MeV 
(similar to the first generation fermion mass) and M ~ 10 9 GeV, we get this eigenvalue to 
be of order 10~ 6 — 2.5 x 1CT 5 eV. Its square is therefore in the range where the VO solution 
to the solar neutrino puzzle can be applied. Also for m ~ a, the mixing angle between the 
zero eigenvalue mode and this mode is maximal. 

In the diagonalization, we have ignored the radiative corrections that allow us to ex- 
trapolate to the weak scale the above mass matrix which is valid at the string scale. Going 
to the weak scale, M z = mass of the Z° meson, there are two kinds of contributions that 
dominate the radiative corrections: one arising from the top quark coupling to the standard 
model Higgs doublet in the effective LHLH operator induced by the seesaw mechanism 
T6| , and a second one which can arise from self couplings of the Higgs fields (for the non- 



super symmetric version of the model). The top quark contribution replaces the parameter 
m m the matrix in Eq. 0bym(l + -^ln(M*/M z )) in the off-diagonal terms, whereas 
the self-scalar coupling contributes only to the m 2 term. Thus the m? term in Eq. |13| is 

replaced by m 2 (l + ^£ln(M*/M z ) + j^l^M* / M z )^j ■ (This expression will have more 
terms involving extra scalar self couplings if there is more than one Higgs doublet in the low 
energy theory). Let us denote the rest of the radiative corrections by e in the m 2 term and e' 
in the off-diagonal ones. Such radiative corrections could come from one loop contributions 
at the scale M itself. The magnitude of the radiative corrections also could increase if the 
low energy theory below M* has more than one Higgs doublet, as noted. It is therefore not 
implausible to assume that the radiative corrections are significant. With redefinition of the 
off-diagonal m terms, M. takes the same form as in Eq. [U| except the m 2 term is replaced 
with m 2 + A, with A = (m 2 /16tt 2 )(24/i 2 A + 4A 2 + e - e')ln(M*/M z ). The characteristic 
equation for the m n becomes 

m 2 + A a 2 7c(m n - A/M) vrm„ 

cot = m„. 15 

M Mfx V 

It is now easy to see that if A << (a 2 +m 2 ), the zero eigenvalue is replaced by ~ M(S+rn?) • 
Clearly, we want A closer to m 2 ,a 2 to obtain our desired parameters. For this purpose, 
we choose parameters A and the other radiative corrections (i.e., e, e') appropriately, so 
that the two lowest eigenvalues are of almost equal magnitude, and the parameters of the 
next section can be reproduced. This situation can be realized more easily if the theory 
is nonsupersymmetric all the way to the string scale. If there is supersymmetry all the 
way down to the TeV scale, then the scalar self coupling contributions "kick" in only below 
Msusy and one has to stretch parameters and perhaps require a two-Higgs doublet model 
below the SUSY scale to realize the parameters used below. 

Let us now turn to the determination of the mixing angles. For this we need the explicit 
form of the eigenvector \l/ n for the n-th mass eigenstate: 



7 



/ 1 \ 



b'l 



hi 



where 



a 
m 



\ ■ l 



(16) 



(17) 



and 



62 = rrm> (18) 



= 3% (19) 



rrir. 



The normalization of the state is given by 
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The mixing of the z/ e with the bulk modes is essentially given by the 1/N n , which for the 
n-th eigen mode (with m n ~ nfi) is~ For ma/M ~ 10 -5 eV and fi ~ 10 -3 eV, this 

mixing is of order of a percent and is therefore in the interesting range for a possible MSW 
transition of the solar neutrinos. 



V. SOLAR NEUTRINO DATA FIT BY A COMBINATION OF VACUUM AND 

MSW OSCILLATIONS 

In this section, we discuss the question of how to understand the solar neutrino data in 
these models, while at the same time explaining the atmospheric as well as the LSND data. 
There have been several recent papers that have addressed the issue of explaining observed 
oscillation data in models with large extra dimensions [|TX)| , pr7| , p~8| . p~5] , p~9| ^23[1 . In particular, 



in Ref. |L5[ it has been shown that while the overall features of the solar and atmospheric 
data can be accommodated in minimal versions of these models with three bulk neutrinos, 
it is not possible to explain simultaneously the LSND observation for the — v e oscillation 
probability, and one must incorporate new physics in the brane. 

In the models presented here, the — v T mass difference responsible for atmospheric 
oscillation data is generated via the radiative corrections in the TeV scale models, and the 
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seesaw mechanism in the local B-L models. Since we arrange the models so that the mass 
of the v^ T pair is in the eV range, this provides a way to accomodate the LSND results. Let 
us therefore focus on the solar data. We will present our discussions using the parameters 
of the TeV scale model. The discussion also applies to the local B-L models, with only the 
labels of the parameters changed. 

A first glance at the values of the parameters of the model such as m in Eq.(3) and 
R" 1 ~ 10~ 3 eV suggests that perhaps one should seek a solution of the solar neutrino data 



in these models using the small-angle MSW mechanism JT7j. However, the present Super- 
Kamiokande recoil energy distribution seems to disfavor such an interpretation, although 
any definitive conclusion should perhaps wait till more data accumulates. In any case if the 
present trend of the data near the higher energy region of the solar neutrino spectrum from 
Super-Kamiokande persists, it is likely to disfavor the small-angle MSW solution and tend 
to favor a vacuum oscillation. However, the latter does not give correct rates for the three 
types of solar neutrino experiments. Note that for sterile neutrinos the large-mixing-angle 
MSW solution does not work. 



As discussed in [I0~|, our solution to the solar neutrino data consists of two components: 



one involving the vacuum oscillation of u e to v B and the second one involving the MSW 
transition of the higher energy z/ e 's to higher KK modes of the vb- The vacuum oscillation 
part is straightforward, and in order to get a better fit we have to adjust the Am 2 _ ■ On 
the other hand, to discuss the MSW effect for the case of bulk neutrinos, we need to include 
the effect of solar matter on the infinite dimensional neutrino mass matrix. 

To have a physical understanding of our strategy, note that the simplest way to reconcile 
the rates for the three classes of solar experiments is to "kill" the 7 Be neutrinos, reduce the 
8 B neutrinos by half and leave the pp neutrinos alone. To achieve this using pure vacuum 
oscillation, one may put a node of the survival probability function P ee around 0.86 MeV. 
However, for an arbitrary node number, the oscillatory behaviour of P ee before and after 
0.86 MeV cannot in general satisfy the other two requirements mentioned above; specifically, 
if there are more nodes prior to 0.8 MeV, the Gallium pp neutrinos get suppressed. If one 
uses the first node to "kill 7 Be", then for 8 B neutrino energies the P ee is close to one and 
not half as would be desirable. The strategy generally employed is not to have a node at 
the precise 7 Be energy but rather somewhat away so that it reduces 7 Be to a value above 
zero. This requires that one must reduce the 8 B neutrinos by much more than 50%, so one 
can fit Chlorine data. The water data then requires an additional contribution, which, in 
the case of active vacuum oscillation (VO), is provided by the neutral current cross section, 
amounting to about 16% of the charged current one. Thus in a pure two-neutrino oscillation 
picture, VO comes close to working for oscillation to active neutrinos but certainly does not 
work for active to sterile oscillation. It is here that the large extra dimensions come to the 
rescue. 

In our model, both vacuum oscillations and MSW oscillations are important. This is 
because the lowest mass pair of neutrinos is split by a very small mass difference, whereas the 
KK states have to be separated by > 10~ 3 eV because of the limits from gravity experiments. 
We can then use the first node of P ee to suppress the 7 Be. Going up in energy toward 8 B 
neutrinos, the survival probability, which in the VO case would have risen to very near one, 
is now suppressed by the small-angle MSW transitions to the different KK excitations of 
the bulk neutrino. This is the essence of our new way to fit the solar neutrino observations 
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In order to compute the effect of solar matter on neutrinos produced deep in the sun, 
we begin with how neutrinos propagate. Eigenvectors of the neutrino mass matrix, At, 
evolve according to e l (p x ~ Et )^ where for a state of mass m, px — Et m E{x — t) — tm 2 /2E. 
Neutrino oscillations happen because different mass states interfere with each other. When 
the neutrino of a particular energy is detected at a particular point, E(x — t) is independent 
of which mass state contributes. The common phase factor ^ %E ^ x ~ t ) can be factored out of 
all states, because we are concerned only with relative phases when considering interference 
terms. Each mass state evolves separately in a vacuum according to 

i d 4 = %S. (21) 
dt IE K ' 

The electron neutrino is an eigenstate of neutrino interactions with matter, but is not an 
eigenstate of the mass matrix, At. In the basis of eigenstates of neutrino interaction with 
matter, an arbitrary state evolves through a vacuum according to 

.da M ] M _ 

i— = — - a. 22 

dt 2E v ' 

In matter the squared mass matrix, At^At, is replaced by AA^AA + 2EHi, where Hi is 
p e = \f2Gp{n e — 0.5n n ) when acting on v e and is zero on sterile neutrinos. The survival 
probability, P ee , is the probability that an initially pure v e is still a v e after the neutrino state 
propagates from its origin in the sun to its detection on the earth. For the TeV scenario, 
P ee depends on E and the parameters used in Air e v of Eq. ||: P ee = P^ e eV (E,5 ee ,m, p). 
But since the propagation depends only on Aft Atj2E + Hi, the survival probability must 
be unchanged whenever M. / \J~E~ is unchanged. The amount of computation required could 
therefore be greatly reduced by use of the scaling rule 

F£ V (E, 5 ee , m, p) = P™ (l, ^, JL^ . (23) 

Similarly, there is a scaling rule for P ee = Pf~ L {E, M, m, A, a, p) of the local B-L scenario: 

P e B ~ L (E, M, m, A, a, p) = Pg~ L ( 1, 1, -^=, -A=, -^=, . (24) 

To discuss the MSW effect for an infinite component system, we first diagonalize the 
matrix Av Ai + 2EH\ for both models. We give the results for the TeV scale model first, 
and in a subsequent paragraph present the result for the local B-L case. 



Te V scenario and matter effect: 



When Ai = Ai.TeVi to express the eigenvectors and eigenvalues of Ai^ A\ + 2EH 1: define 
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(25) 



Wk = 1 in vacuum. The characteristic equation becomes 
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(26) 



The eigenvectors are as in Eq. || except the coefficients of z/ s and u' B _ acquire an additional 
factor l/w n , the m 2 , is replaced by m 2 ,, and the normalization becomes 



Nl 



+ f 2 2n2 , (m n -wAe) 2 ' 

2 I m 2 



( X J^) (m n - w n 5 ee ) 
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(27) 



Note that in the presence of a dense medium, very crudely speaking, the p e term in Eq. ^ 
will dominate over the rest of the terms, and when m 2 ~ k 2 p 2 for any of the KK levels, that 
particular coefficient in Eq. H dominates, and the MSW resonant condition is satisfied. 



Matter effect in the local B-L case: 



Following the same procedure, we get for the local B-L case the following eigenvalue 
equation in the presence of matter effects (where we will ignore the radiative corrections 
since they do not affect the results materially): 
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(28) 



If we denote the eigenvector of the matter-affected mass matrix as 



we find that 
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* = m 

Here again we see that when the density term dominates rfi|, there can be an MSW transition 
to the nth level when m 2 k ~ n 2 /i 2 . 

To carry out the fit, we studied the time evolution of the v e state using a program that 



evolved from one supplied by W. Haxton [24]. The program was updated to use the solar 



model of BP98 |25|] and modified to do all neutrino transport within the sun numerically. For 
example, no adiabatic approximation was used. Changes were also necessary for oscillations 
into sterile neutrinos and to generalize beyond the two-neutrino model. Up to 16 neutrinos 
were allowed, but no more than 14 contribute for the solutions we considered. While we 
explored the parameter space using BP98, the more recent BP2000 [EE] solar model gave 



almost identical results where the two models were compared. 

For comparison with experimental results, tables of detector sensitivity for the Chlo- 
rine and Gallium experiments were taken from Bahcall's web site [E3[]. Neutrino spectra 



from the various solar reactions were taken from the same site, except for the B spectrum, 
which comes from a more recent determination ||27|| . The Super-Kamiokande detector sen- 



sitivity was modeled using |28f| , where the percent resolution in the signal from Cherenkov 
light, averaged over the detector for various total electron energies, is provided. To within 
the number of digits accuracy given, the fractional resolution is (.443 + .0038E)/\^E, 
where E is the total electron energy. Combined with knowledge of the relation between 
amount of Cherenkov light and the true electron energy, this gives an energy resolution of 
ge = E ( E ~ E ^^ u 3y)_oo38E) ^ rjn^ e reS p 0nse Q f Super-K to neutrinos of energy E u is given by 



smearing the differential cross section given by 't Hooft (including radiative corrections 



accounted for by modifications of qa and gy [30]) with a Gaussian resolution function of 
standard deviation erg, and restricting the integration to measured energy within an energy 
bin, or within the range of the total flux measurement (5.0-20 MeV). 

Calculations of electron neutrino survival probability, averaged over the response of de- 
tectors, were compared with measurements. While theoretical uncertainties in the solar 
model and detector response were included in the computation of x 2 as described in Ref. 



[51 1, the measurement results given here include only experimental statistical and system- 
atic errors added in quadrature. The Chlorine survival probability, from Homestake |32|, is 



0.332 ± 0.030. Gallium results |5| for SAGE, GALLEX and GNO were combined to give 
a survival probability of 0.579 ± 0.039. The 5.0 — 20 MeV 1258 day Super-K experimental 
survival probability [|34| is 0.451 ± 0.016. The best fits were with R « 58/im, mR around 



0.0094, and 5 ee ~ 0.84 x 10~ 7 eV, corresponding to 5m 2 ~ 0.53 x 10~ n eV 2 . These param- 
eters give average survival probabilities for Chlorine, Gallium, and water of 0.383, 0.533, 
and 0.450, respectively. They give a v e survival probability whose energy dependence is 
shown in Fig. |l|. For two-neutrino oscillations, the coupling between v e and the higher mass 
neutrino eigenstate is given by sin 2 26, whereas here the coupling between v e and the first 
KK excitation replaces sin 2 26 by 4m 2 R 2 = 0.00035. 
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FIGURES 



1 r- 

0.9 - 




FIG. 1. Energy dependence of the v e survival probability when R = 58/um, mR = 0.0094, 
5 ee = 0.84 x 10~ 7 eV. The dot-dashed part of the curve assumes the radial dependence in the Sun 
for neutrinos from the pp reaction, the solid part assumes 15 radial dependence, and the dashed 
part assumes 8 B radial dependence. 

Vacuum oscillations between the lowest two mass eigenstates nearly eliminate electron 
neutrinos with energies of 0.63 MeV/(2n+l) for n = 0, 1, 2, ... . Thus Fig. [I] shows nearly 
zero u e survival near 0.63 MeV, partly eliminating the 7 Be contribution at 0.862 MeV, and 
giving a dip at the lowest neutrino energy. Increasing 5 ee moves the low energy dip to the 
right into Gallium's most sensitive pp energy range, making the fit worse. Decreasing S ee 
increases Gallium, but hurts the Chlorine fit by moving the higher energy vacuum oscillation 
dip further to the left of the 7 Be peak. Note that the pattern of two eigenstates very close in 
mass persists for the Kaluza-Klein excitations as well. These MSW resonances start causing 
the third and fourth eigenstates to be significantly occupied above ~ 0.8 MeV, the fifth 
and sixth eigenstates above ~ 3.7 MeV, the 7'th and 8'th above ~ 8.6 MeV, and the 9'th 
and 10'th above ~ 15.2 MeV. Fig. [I] shows dips in survival probability just above these 
energy thresholds. The typical values of the survival probability within the 8 B region (~ 6 
to ~ 14 MeV) are quite sensitive to the value of mR. As can be seen from Eq. [|, higher 
mR increases 1/N n « m/m n « mR/n for various n, and thereby increases u e coupling to 
higher mass eigenstates, strengthens MSW resonances, and lowers v e survival probability. 

The expected energy dependence of the v e survival probability is compared with Super-K 
data |34| in Fig. 0. The uncertainties are statistical only. The parameters used in making 
Fig. ^| were chosen to provide a good fit (x 2 = 3.4) to the total rates only; they were not 
adjusted to fit this spectrum. But combining spectrum data with rates using the method 
described in Ref. gives x 2 — 14.0 for the spectrum predicted from the fit to total rates. 
With 18 degrees of freedom, the probability of x 2 > 14.0 is 73%. If instead the fit were to 
an undistorted energy spectrum, the x 2 would be 19.0. 
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5 6 7 8 9 10 11 12 13 14 

Recoil Electron Energy (MeV) 

FIG. 2. Super-Kamiokande energy spectrum: measured results based on 1258 days (error 
bars) and predicted (curve) for the same parameters as in Fig. 1. The curve is not a fit to these 
data. 

One may also seek fits with 5 ee constrained to be very small, thereby eliminating vacuum 
oscillations. We found x 2 — 4.4 for the best such fit to the total rates only. The same param- 
eters then used with the Super-K spectrum gave x 2 — 18.7, corresponding to a probability 
of 41%. 

The seasonal effect was computed for a few points on the earth's orbit. If r is the distance 
between the earth and the sun, 

— = l + ecos(6 -6 Q ), (31) 

where r is one astronomical unit, e = 0.0167 is the orbital eccentricity, and 8— 6 « l-nit— £ ), 
with t in years and to = January 2, 4h 52m. Table | shows very small seasonal variation. 
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TABLES 

TABLE I. Predicted seasonal variations in v e fluxes, excluding the 1/r 2 variation. The model 
assumed fi = 0.32 x 1(T 2 eV, m = 0.34 x 10" 4 eV, and 5 ee = 0.78 x 10~ 7 eV. 



9 — 9q in eqn 31 


Chlorine 


Gallium 


Water 


(January 2) 


0.3787 


0.5144 


0.4635 


±vr/2 


0.3762 


0.5121 


0.4633 


7r (July 4) 


0.3747 


0.5082 


0.4631 
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VI. CONSEQUENCES AND COMMENTS 



While the seasonal effect is too small to be observable soon, and day-night effects would 
also be difficult to see, the characteristic energy spectrum of Fig. [I] should be observable 
by SNO, the first results from which will be available shortly. Whereas the characteristic 
dips in energy are nearly washed out in the neutrino-electron scattering results of Super- 
Kamiokande shown in Fig. [|, they should be seen clearly with the far better resolution of 
SNO's charged-current interaction. If that confirmation of the fit presented here eventuates, 
then the resulting 0.06 mm extra dimension size should be directly detectable by gravity 
experiments in the not too distant future, since the present best limit on such effects is 
less than a factor of four from that value. 

In contrast, the mass of the electron neutrino, which consists mainly of eigenstates of 
mass 3 x 10~ 5 eV, is undetectable directly by such means as tritium endpoint experiments. 
Whle neutrinoless double beta decay measures an effective neutrino mass, the additional 
contributions from and v T are probably so small as to make that process unlikely to be 
observed, although there are conjectured mechanisms other than neutrino mass to produce 
neutrinoless double beta decay. 

If the relatively large dimension size is confirmed by SNO, this will raise issues about 
cosmological and supernova limits from the effects of high Kaluza-Klein states of both sterile 
neutrinos |18|JT9| , |37[| and gravitons |38] . These constaints are much stronger than those from 



new contributions to low energy weak processes |5I| , but are somewhat suspect because the 
two regimes, the early universe and the supernova core, are very complex and not yet fully 
understood. Also, for sterile neutrino limits, the phenomenology presented here is aided 
by there being a single Kaluza-Klein tower based on an exceedingly small mass, the VO 
Am 2 is an order of magnitude smaller than usual, and for MSW the equivalent sin 2 2# is 
more than an order of magnitude smaller than for standard fits. For the global B-L model, 
the universe re-heat temperature could be very low, since anything above 0.7 MeV has 
cosmological validity, reducing production of high KK states. The high string scale of the 
local B-L model appears to avoid all these constraints, however. This suggests that more 
complete investigation of the constraints may provide a means of choosing between these 
quite different models [^TJ both of which can provide this new way to give a simultaneous 
fit to the solar, atmospheric, and LSND data, possibly giving the first evidence for an extra 
large dimension. 

The work of R.N.M. is supported by a grant from the National Science Foundation No. 
PHY-9802551. The work of D.O.C. and S. J.Y. is supported by a grant from the Department 
Of Energy No. DE-FG03-91ER40618. 
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Appendix A: 

Model for the Majorana mass matrix for the brane neutrinos 



In this appendix, we seek a possible theoretical origin of the neutrino mass pattern used 
in the TeV scale model. We will keep the standard model gauge group and attempt to 
extend the higgs sector in such a way that one generates the neutrino masses at loop levels. 
The basic idea will be to consider L e + L^ — L T symmetry for neutrinos, with L e — 1 for vb- 
On general symmetry grounds, the allowed mass matrix for (u e , u^, v T ) can be written as 

/ S eT \ 

M v = m . (32) 
\5 er m / 

The remaining terms are assumed to arise after we turn on the symmetry breaking so that 
those elements are small. 

For an explicit realization, we augment the standard model by the singlet charged Higgs 
field h ++ , which is blind with respect to lepton number, h' e , and which carries L e = —2, and 
SU(2)l triplet fields, A ejMjT , with Y = 2, which carry two negative units of lepton numbers 
L e ,fi,r, respectively. The Lagrangian involving these fields consists of two parts: Co and C\. 
Co is invariant under (L e + L M — L T ) number and is given by: 

C = ferh ++ e~ R r R + f^ T h ++ ^ R T R + f' ee h ++ 'e R e R 

+g e L e L e A e + g^L^L^A^ + g T L T L T A T + h.c. (33) 

We assume the Higgs potential to contain the L e + — L T invariant terms 

V = £ M h ++ (A 2 a ). (34) 

a=e,fj,,r 

The soft breaking of the symmetry of L e + L^ — L T is achieved by the following terms in the 
Lagrangian: 

A = h ++ ( E M *A 2 + Mo^A T + M 0r A e A r ) + ii 2 h ++ ^h ++ ' + h.c. (35) 

i=e,/j.,T 

With these couplings, the neutrino Majorana masses arise from two-loop effects (similar to 
the mechanism of ref. |4^|), and we have the L e + L^ — L T violating entries Sij ~ c m e m e .. 
Using 5 ee ~ 10~ 8 eV, then 5 TT ~ 10~ 2 , as would be required to understand the atmospheric 
neutrino data. In Fig. we give a typical two loop diagram that leads to neutrino masses. 
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FIG. 3. Typical two loop diagram for neutrino mass in the TeV string scale model 

Appendix B: 
Neutrino mixings in the local B-L model 



In this section, we discuss the details of neutrino mixing in the local B-L model described 
in the text (section 3). In the text we considered only the v e — v B sector. However, in order 
to explain LSND results, we have to generate the v e — mixing and make sure that its 
presence does not affect the considerations in the text. 

For this purpose, we consider a model based on the gauge group SU{2)l x SU{2)r x 
U{1)b-l with quark, Q = (u,d) and lepton L = (z/ e ,e) doublets assigned as usual in 
a left-right symmetric manner and Higgs fields belonging to bidoublet field <j)(2, 2, 0) and 
B — L = 1 doublets Xl,r- In addition, we require the theory to respect a softly broken 
global V = L e — + L T symmetry. The t/(l)z/ invariant part of the Lagrangian can be 



written as: 



C 



L {ll>a,RXR){ll>b,RXR) + {R-^ L ) _|_ I, [(^e,RXR)(j>e,RXR) + (R-^L)]8 

+ Jdy v B l(x, y)d 5 v B R(x, y) + h.c, 



(36) 



where the form of the matrix is determined by the U(1)l> symmetry. After symmetry 
breaking, i.e. (5) ~ (xr) 7^ 0, Eq. ^6| leads to a right-handed Majorana neutrino mass 
matrix of the form: 



M, 



R 



( M ee M e(i 



M, 



efj. 



\ M , 



M, 







The complete mass matrix for 



fl,T 



(37) 



> e , u T , N e , Np, N T , v B ) is given by 
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(38) 



Let us now consider a hierarchical form for the Dirac masses i.e. vn\ <C m 2 <C m 3 . Looking 
at the "non-bulk" part of this mass matrix, one can see that if we choose M efl ~ M ee ~ 
and m 2 ~ m 3 /10 ~ 100 MeV, then one gets the right value for the Am\ TMOS and the v e — 
mixing angle needed for understanding the LSND results. Furthermore, one can decouple 
all the heavy neutrinos as well as the z/ M T from the spectrum and obtain the same equation 



as in 13 and all the considerations in the B-L model described in the text go through. 



Appendix C: 

Naturalness of ultra light sterile neutrinos in brane-bulk models 



The bulk neutrino "self mass" terms are constrained by the geometry of the bulk and 
could therefore under certain circumstances be zero. If that happens, the only mass of 
the KK states of the Vb will arise from the kinetic energy terms such as vb^ i Qi^b^ where 
/ = 5,6... and will be given by n/R for R the radius of the extra dimensions. In such a 
situation, an ultralight vb,kk arises naturally. 

The key to naturalness of the ultralight bulk neutrino is the geometry that forbids both 
Dirac and Majorana mass terms. Let us give a few examples. In five dimensions, if we impose 
the Z 2 orbifold symmetry (y — > —y), then it follows that the Dirac mass vanishes. Now if 
we impose lepton number symmetry in the brane, the Majorana mass vanishes, leaving us 
with no mass term for the bulk neutrino in 5-dimensions. 

Another interesting example is the 10- dimensional bulk, where the bulk neutrino is a 
16-component spinor which, when reduced to four dimensions, leads to eight 2-component 
spinors. The interesting point is that for a 16- dimensional spinor, one cannot write a Dirac 
or Majorana mass term consistent with 10-dimensional Lorentz invariance. In this case, 
there is no need for assuming lepton number symmetry to get an ultralight sterile neutrino. 
A similar situation is also expected in six dimensions, if we choose the bulk neutrino to be 
a 4-component complex chiral spinor. 
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